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By considering both the average structural parameters obtained from Rietveld refinement of 
neutron powder diffraction data, and the local structural parameters obtained from the atomic pair 
distribution function, we have tested the recent hypothesis of Perry et al. [Perry et al., Phys. Rev. B 
65, 14450f (2002)] that doping in La2-a;Sri:Cu04 system occurs as localized defects of predominantly 
Cud322_r2-Opz character associated with the Sr dopants accompanied by a local destruction of 
the Jahn- Teller distortion. While the structural parameters behave qualitatively according to the 
prediction of this model, a quantitative analysis indicates that doped holes predominantly appear 
in the planar Cud^2 _y2-Opx,y band as is normally assumed. However, a small amount of the 
doped charge does enter the C\xd^^2 _,.-2-Opz orbitals and this should be taken into account when 
theoretical phase diagrams are compared to experiment. We present a calibration curve, p = 
2;(1.00(f ) — 0.45(7)a;), for the planar charge doping, p, vs strontium content, x, for the La2_iSr^Cu04 
system. 

PACS numbers: 74.81.-g,74.72.Dn,74.72.-h,61.12.-q 



I. INTRODUCTION 

Cuprate high temperature superconductors are doped 
Mott insulators. The novel superconductivity appears at 
doping levels just beyond the insulator-metal (IM) tran- 
sition. The insulating behavior of the undoped endmem- 
ber is understood to be due to electron correlation ef- 
fects in the half-filled planar Cxidx^-y'^-Opa (o: = x,y) 
anti-bonding band. The phase diagram of the cuprates 
can be interpreted in terms of holes doped into this pla- 
nar band. For example, in the system La2-a;Srj;Cu04 
it is assumed that one hole enters this band per stron- 
tium atom. This paradigm for the doping has hardly 
been questioned and a multitude of papers exist implic- 
itly assuming this behavior. The La2-a:Sr3:Cu04 system 
is archetypal since it is a single layer system which can be 
straightforwardly doped over a wide range. The result- 
ing phase diagram is thought to exhibit features that are 
universal to the cuprates. However, as a result of recent 
ah initio electronic band structure calculations. Perry, 
Tahir-Kheli and Goddard (PTG)^ have suggested a new 
model for the doping in the La2-a;Sra;Cu04 system. In 
this case the doped holes reside in localized states on 
the CuOg octahedra situated next to the dopant stron- 
tium ions. Furthermore, the doped charge resides prin- 
cipally in the G\xd^z^-r^-Opz orbitals (i.e., the out-of- 
plane bonds). Coincidentally a structural distortion oc- 
curs such that the local Jahn- Teller (JT) distortion, that 
results in the long Cu-0 apical bond, is destroyed. A flat 
impurity band forms in the gap. The IM transition then 
occurs on increased doping in a manner similar to doped 
semiconductors by percolation of the doped impurities. 
These calculations were originally motivated by XAFS 
measurements that observed a distorted environment of 
the octahedra in the vicinity of dopant ions'^. This new 
way of understanding the doping, if it is right, clearly 
will result in a paradigm shift in our understanding of 



cuprate physics. It is thus of the greatest importance 
to test the hypothesis experimentally. A rather direct 
probe of this doping mechanism is the structure because 
it involves the local destruction of JT elongated bonds. 
The average Cu-0 bond length along z will thus shorten 
with doping and the width of the bond-length distribu- 
tion for this bond will increase reflecting the increased 
disorder. Despite detailed structural studies of this sys- 
tem^ the temperature and doping dependence of these 
parameters has not been reported. In this paper we in- 
vestigate whether there is evidence supporting this new 
doping paradigm in both the average crystal structure, 
and the local structure as measured using the atomic 
pair distribution function (PDF) technique^. We reex- 
amine the extensive earlier structural data of Radaelli to 
extract the desired parameters. We have also analyzed 
new neutron powder diffraction data using Rietveld re- 
finement and PDF refinement. We find that the data 
qualitatively agree with the predictions of the Perry et 
ali: However, a more quantitative analysis, including 
evaluating the bond valences as a function of doping of 
the in-plane and apical Cu-0 bonds, indicates that doped 
charge is predominantly residing in the planar bonds sug- 
gesting the existing paradigm for doping is valid. 



II. EXPERIMENTAL 

This study reexamines the extensive published data of 
Radaelli et alA to extract as a function of doping the in- 
teresting parameters of the plane copper to apical oxygen 
bond length, rcu-02 and the short lanthanum/strontium 
to apical oxygen bond length, ?'La/Sr-02- These param- 
eters were determined from the reported fractional co- 
ordinates and lattice parameters. Also important is the 
doping dependence of the copper to in-plane oxygen bond 
length, Tcu-oii which was presented previously^ but is 
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FIG. 1: Representative full profile PDF fit of standard LTO 
model to 10 K data for undoped sample. Experimental profile 
shown as open circles, model profile as solid line. Difference 
curve is shown below the fit as solid line. Dashed lines denote 
experimental uncertainties at the 2a level. PDFs studied here 
utilize diffraction information up to Qmax =35 A ~^ . 



reexamined here. To check the PTG prediction we also 
need the anisotropic displacement parameter along the 
z-axis of the apical oxygen, [733(02). This was not pub- 
lished in^ so we have used recently collected neutron pow- 
der diffraction data of our own from a less extensive set 
of samples. Samples were synthesized using solid state 
methods and loose powders of 10 g sealed in vanadium 
cylinders were measured at 10 K a the GEM diffractome- 
ter at the ISIS neutron source. Details of sample synthe- 
sis, characterization and measurement are reported else- 
where^. Rietveld refinements were carried out on the 
data using the program GSAS-. The data were also cor- 
rected for experimental effects and normalized to obtain 
the atomic pair distribution function^ (PDF) using the 
program PDFgetN?^. The PDF is the Fourier transform 
of the normalized corrected powder diffraction data. It 
utilizes both Bragg and diffuse scattering and contains in- 
formation about the local structured. Structural mod- 
els were refined to the PDFs using the program PDF- 
FIT—. To check for consistency, the same parameters 
were varied in the PDF and Rietveld refinements. An 
example of the measured PDF from the undoped sample 
is shown in Figuren]with the best fit model-PDF plotted 
on top. The PDF allows us to explore models that con- 
tain non-periodic defects such as those proposed by PTG. 
The PDF was calculated from a supercell model with the 
doping induced defects discussed by PTG. This was then 
used as simulated data and fit using the usual long-range 
orthorhombic (LTO)^ structural model. In this way it 
was possible to estimate the size of the enlarged displace- 
ment factors that would result from the presence of the 
doping induced disorder described by PTG. In Perry et 
al} , for simplicity, no octahedral tilts were considered. In 
the present case we must compare resulting models with 
real data and so the doping induced defects described ini 
were superimposed on the background of octahedral tilts 
observed in the undoped endmember. 




0.1 0.2 0.3 0.4 

strontium content, x 



FIG. 2: Evolution of rcu~02 and ri^a/Sr-02 (a), and rcu-o\ 
(b) average bond lengths with hole doping at 10 K: Rietveld 
result from^ (solid circles) and PDF result (solid squares). 
Dotted lines: slope predicted from simple electrostatics con- 
siderations. Dashed lines: slope prediction based on PTG 
model. Arrows (from left to right): IM transition, structural 
phase transition, and disappearance of superconductivity. 



III. 



RESULTS AND DISCUSSION 



The qualitative predictions for the evolution of the av- 
erage structure with doping in the scenario of Perry et 
al.i- are: (1) decrease in the average rcu-02 due to 
the local destruction of Jahn- Teller distortions at Cu'^"'" 
sites; (2) increase in the [733(02) with doping because of 
the coexistence of Jahn- Teller distorted and non-Jahn- 
Teller distorted octahedra; (3) increase in the average 
'i'La/Sr-0{2) distancc. All these effects are seen qualita- 
tively in the data as shown in Figs. |21 and El 

However, similar effects can be expected even if the 
doping is taking place in the conventional way: homoge- 
neously into the planar bonds. For example, if the Cu02 
planes are becoming more positively charged with doping 
the negatively charged apical oxygen would be predicted 
to come closer to the planes due to simple coulomb at- 
traction. Also, the U33 displacement factor of the apical 
oxygen atom is expected to increase with doping due to 
the random doping of larger Sr^"'' ions. Inhomogeneous 
doping in the form of stripes or checkerboards ^i^iSi would 
also result in enlarged thermal factors for in-plane and 
apical motions of Ol and 02. However, inhomogeneous 
doping would not affect average bond-length obtained 
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from Rietveld refinement. It is therefore important to be 
more quantitative to distinguish these different possibil- 
ities. 

The PTG calculations^ were carried out at special ra- 
tional doping fractions of 1/8, 1/4 and 1/2 and suggest 
the appearance of a single distorted CuOe octahedron as- 
sociated with each doped strontium. The rcu-02 bond 
closest to the strontium is shortened by Arcu-02 = 
—0.24 A and that farthest from the Sr on the same octa- 
hedron by Arcu-02 — —0.10 A, the Sr-02 distance in- 
creases, Ari^g./Sr-02 = +0.11 A and the in-plane tqu-oi 
distances do not change. In the average structure these 
defects are not seen explicitly; however, the distortions 
will be apparent as a properly weighted change in the 
average bond-length and an increase in the respective 
displacement parameters. If we assume that these de- 
fects appear as each strontium is doped and are not a 
special feature of the rational doping fractions studied, 
we get the following relations for the strontium doping, 
X, dependence of the average bond lengths: 

Arcu-02{x) = -0.17 X, ArLa/Sr^02{x) = 0.11 X, (1) 

where x is the doping level. These are shown as the 
dashed hues in Fig. |21 The evolution of rLa/Sr~02 is 
in quantitative agreement. The case of rcu-02 is more 
complicated. rcu-02 decreases much more slowly than 
the prediction initially on doping. However, beyond a 
doping level of 0.17 the slope of rcu-02 vs x increases. 
In this higher doping region the slope of the PTG line 
agrees rather well with the data. There is no specific pre- 
diction from the PTG calculations for the x-dependence 
of rcu-oi, although it is expected to be small since the 
length of rcu-oi in the defect octahedra does not change 
and no doping dependent change in lattice parameters 
has been assumed in their calculations. This reflects the 
fact that, in the PTG picture the doped charge is almost 
exclusively located in the apical d^2_^2 orbitals of the 
copper. Therefore, we have drawn a flat dashed line as 
the PTG prediction for this parameter in Fig. [2Ib). As 
is apparent, the data for rcu-oi slope downward with 
a slope of Arcu-oi ~ —0.105. This is comparable to, 
though slightly less than, the observed Arcu-02 ~ —0.18 
of the apical bond. This observation will be important 
later. 

The presence of the anti-Jahn- Teller distorted doped 
CuOg octahedra in the PTG calculations will manifest 
itself in an increased atomic displacement factor for 02, 
especially along z. The most sensitive parameter will 
be the width of the Cu-0 pair distribution along the z- 
axis. This cannot be measured directly from the PDF 
because of peak overlap but can be found from the sum 
of U33(Cu)-|-U33(02). The mean-square width is plot- 
ted in Fig. |31 and does increase with increasing doping 
at constant temperature, a sign that static disorder as- 
sociated with this bond is increasing with doping. To 
investigate more quantitatively whether this increase is 
consistent with the PTG predictions we calculated the 
PDF of a model structure containing the PTG defects. 
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FIG. 3: Sum of displacement parameters U33 for Cu and 02 
as a function of strontium content at 10 K, from PDF refine- 
ments. Solid line is a guide to the eye. Dashed line denotes 
expected increase in the PDF displacement parameters, es- 
timated from PDF simulation based on PTG model for 1/8 
doping (see text). 



This simulated PDF was then fit using PDFFIT using 
the standard LTO model. This gives us a measure of 
how much {733(02) would increase in a standard LTO- 
structure refinement due to the presence of the PTG de- 
fects. The PDF was simulated for the x = 0.125 com- 
position using a 8x unit cell in the following way. The 
reference structure used was that of the undoped mate- 
rial with atomic positions taken from^ and displacement 
factors taken from our PDF result for the undoped sam- 
ple. One La in the supercell was then replaced with Sr 
and the atoms on the CuOg octahedron adjacent to the 
Sr were displaced by the A-values from the PTG predic- 
tion given above. The displacement parameters in the 
model retained their values in the x = compound. 
This neglects strain due to the defects, and any stripe 
or checkerboard doping inhomogeneities in the Cu02 
planes, that would tend to increase displacement param- 
eters further. The resulting PDFs were refined with the 
standard LTO unit cell. This resulted in an enlarged 
value of U33(Cu)-hU33(02)= 9.1 x 10^^ ^-2 plotted on 
Fig. 121 which overestimates the observed displacement 
factors. 

Next, using PDF data of 1/8 doped La2-a;Sra;Cu04 at 
10 K, explicit modeling of the anti-JT (A-JT) model was 
performed. The outcome is compared to the results of 
the modeling of the standard LTO model. The model- 
ing employed the program PDFFIT^*'. A large supercell 
model was generated that incorporates explicit Sr dop- 
ing, and allows for defining special doping-affected oc- 
tahedral units, as described above in accord with Refi. 
Three distinct variants of the explicit A-JT models were 
attempted. These variants are distinguished as follows: 
(1) A model with symmetric A-JT distortion on affected 
octahedral units (apical Cu-0 distances restricted to have 
the same length). (2) A model with asymmetric A-JT 
distortion (relevant apical distances are allowed to vary 
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TABLE I: Standard LTO model versus A-JT defect models: 
PDFFIT modeling results summarizing relevant distances. 
All lengths are in A. Rw is the weighted PDFFIT agreement 
factor. 





LTO 


A-JT(l) 


A-JT(2) 


A-JT(3) 


Rw 


0.127 


0.122 


0.131 


0.127 


r(Cu-02) 


2.4179(24) 


2.4095(25) 
2.460(19) 


2.4119(22) 
2.450(20) 
2.469(19) 


2.406(3) 
2.49(4) 
2.51(3) 


r(La-02) 


2.346(3) 


2.354(3) 


2.354(3) 


2.361(4) 


r(Sr-02) 


2.346(3) 


2.305(20) 
2.27(3) 


2.305(20) 
2.27(3) 


2.27(4) 
2.20(5) 



independently), and (3) A model with asymmetric A- 
JT distortion this time restricted to have displacements 
along LTO c-axis only. In all three variants atomic pa- 
rameters of Sr were decoupled from those of La, except 
for the displacement factors that were restricted to be 
the same. In the variant (3) Sr y-coordinate was set to 
zero and Sr motion was restricted to be along c-axis only. 
The A-JT models were refined from the same set of ini- 
tial values as that used for modeling the standard LTO 
structure. This equality of the starting values was only 
violated for 02 and Sr fractional y-coordinates in the 
variant (3), where these were set to zero and fixed. 

The results are summarized in Table ^ where the 
weighted PDFFIT agreement factors (Rw)^" are re- 
ported, as well as all relevant distances. Although for all 
the variants of the A-JT model the agreement factor is 
similar to the one for the LTO model, and for some vari- 
ants even somewhat more favorable, it is not sufficient to 
unambiguously favor the distorted model over the stan- 
dard LTO structure. Moreover, the refined distortions 
obtained from such explicit modeling yield relevant dis- 
tance changes in the direction opposite to that originally 
proposed within the A-JT model of Perry et alX (i.e., the 
Cu-02 distance of the affected octahedral unit lengthens, 
while rsr-02 shortens). However, inspection of the dis- 
tances for the A-JT model results in the proper weighted 
average values when compared to corresponding average 
distances obtained by employment of the standard LTO 
model. 

We now consider the expected structural changes that 
would occur based on electrostatic considerations in the 
conventional doping model. First, we note that the 
charge of the Cu02 plane is becoming less negative with 
doping which could result in the apical 0^~ moving closer 
to the plane. A very rough estimate of this behavior is 
possible neglecting the Jahn- Teller effect and assuming 
an ionic picture by using ionic radii of O^^ = 1.35 A, 
Cu2+ = 0.73 A, Cu3+ = 0.54 A, l.a?+ = 1.16 A, and 
Sr2+ = 1.26 Ai^. From this we get Arc„-02 = -0.19 x 
and Ari£i/5,._Q2 = 0.1 x which appear in Fig.|21as dotted 
lines. Not surprisingly, these curves are similar to those 
predicted by PTG that is also a model involving a mix- 



ture of Cu^"*" and Cu^+ sites. Uniform doping of charge 
into the Cu-0 planes would not result in an increase in 
?733(02) on its own; however, the presence of misfitting 
Sr^+ dopant ions would. Any inhomogeneous charge dis- 
tribution in the plane would increase ?733(02)The fact 
that the observed increase in 1/33(02) exists but is smaller 
than needed to explain the PTG defects may argue in fa- 
vor of electrostatic interpretation. There is a suggestion 
from single crystal refinements compared with empirical 
potential calculations that increases in displacement fac- 
tors with doping can be explained in this way without 
invoking PTG-type defects^*. Our results would tend to 
support this. 

A more refined empirical framework for studying the 
distribution of charge between bonds is the bond valence 
model In this theory there is a direct relationship 
established between bond-length, r, and the amount of 
charge in a bond (the bond- valence, s). The bond valence 
is defined as 

s(r)=exp^!:^L_l), (2) 

where B = 0.37 A is a universal constant and tq depends 
on the chemical identity of the ions in question. Using 
this approach it is possible to determine how much of 
the doped charge is going into the planar vs. the api- 
cal bonds from a measurement of their length changes. 
The theory is successful in unstrained compounds and 
has recently been applied to mixed valent crystals and 
crystals containing Jahn- Teller distortionslfe. The nomi- 
nal valence, Vi, of an ion is then determined by summing 
the bond-valences over all the bonds in which it partic- 
ipates, Vi = ^jSij. In the case of mixed valent crys- 
tals some ambiguity exists as to exactly what value to 
use for r{). The qualitative results of our analysis, how- 
ever, do not depend on what value of tq is used. The 
results are shown in Fig. 0] Here we plot how the total 
doped charge (assumed to be given by x) distributes it- 
self between the apical and planar orbitals respectively, 
as determined from the bond-valence sums. The dimen- 
sionless variables, rjp (solid circles) and rja parameterize 
this partition where rjp + rja = ^■ 

Details of these calculations are as follows. First we 
determine the change in bond-valence of a particular 
bond due to doping, 8si{x) = Si{x) — Si(0). We are 
interested in the partition of charge between the planar 
and apical orbitals, therefore we define an orbital- valence 
5Sm — X^i where the sum is over the planar bonds 
for the Gudx2_y2-Opx,y orbital, SSp, and over the api- 
cal bonds for the Gnd^^'^ _r2-0pz orbital, 5Sa- This is 
then expressed in dimensionless parameters by dividing 
by the total excess valence of copper determined by the 
sum over all the orbital- valences, r]m{x) = 5Sm{x)/VQ^. 
What is plotted in Figure^ is xrjp (solid circles) and x'qa 
(solid squares) , that are measures of the amount of charge 
doped into the planar and apical orbitals, respectively. 

As is apparent in Fig. |21 both the apical and in-plane 
bonds shorten with doping indicating a reduction in 
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is not supported by the structural data. To a rather 
good approximation doped charge is going into the pla- 
nar Cvidx2_y2-0px,y orbitals. We do not address here 
whether Cu sites contain different amounts of charge 
and this does not rule out localized doped defects or 
striped or checkerboard patterns of charge if they have 
sufficient planar character. Such defects have been stabi- 
lized in the unrestricted Becke-3-Lee- Yang-Parr density 
functional calculationsii. However, the observation that 
rcu-02 is getting shorter also requires that some doped 
charge is appearing in the apical Cud^^^ _^2-0pz orbital 
and a correction should be made to the canonical p = x 
relationship to account for this, where p is the doped 
charge in the planar orbitals. The partitioning param- 
eter rip gives the correction, and p is therefore given by 
p = xrjp {rjp shown in Fig. ^a)). This augments spectro- 
scopic measurements that show the partition of charge 
between copper and oxygen but struggle to differentiate 
partition between planar and apical bondsiSiiSiSS. 

In Figure Efa) we show rjp on an expanded scale. To 
obtain an analytic form for the p{x) calibration curve we 
have fit curves to the data in Fig.^a). Linear fit proved 
adequate resulting in r/p = — 0.45(7)x -t- 1.00(1). This 
results in the calibration 



FIG. 4: (a) Partitioning parameter vs. Sr content, (b) Esti- 
mated average distribution of doped charge at 10 K based on 
bond valence calculations using Rietveld obtained distances: 
net apical share (solid squares), net planar share (solid cir- 
cles), total doped charge (open circles). Inset: a measure of 
the amount of charge doped into the apical orbitals. Dashed 
lines represent fit to the data (see text for details). 



charge and doped holes going into both planar and api- 
cal bonds. However, one of the key aspects of this theory 
is the non- linear form of the bond valence, s{r). As a 
result, the same amount of charge doped into a long- 
bond shortens it much more than it would a short-bond, 
and vice versa. At x = rcu-oi < fCu-02 refiecting 
the fact that initially the holes that make Cu in the 2+- 
state reside predominantly in the planar orbitals making 
these the canonical half-filled anti-bonding bands, with 
the apical bonds being doubly occupied with electrons 
(i.e., full, no holes). The physics of this is the existence of 
the Jahn- Teller distortion; however, this shows that the 
phenomenological bond- valence model is making the cor- 
rect prediction for the charge segregation initially. Sub- 
sequently, the experimental observation is that both the 
planar and apical bonds shorten at about the same rate 
but because the planar bonds are initially shorter, and 
the fact that there are four of them rather than two, 
mean that the experimental observations demand that 
the doped charge is going predominantly into the planar 
orbitals as is evident in Fig. ^ 

The main result of this analysis is that the simple 
picture of localized Cud^z2 _j,2~0pz doped holes associ- 
ated with Sr sites described by the PTG calculations 



p(a;) = x(1.00(l) - 0.45(7)a;). (3) 

The resulting fit to the data is shown in Fig.^a) and (b) 
and in the inset as a dashed line. Thus, at a strontium 
content of x = 0.125 the doping in the planar orbitals 
is actually p — 0.119, and conversely, the rational 1/8 
filling occurs at x = 0.131. Interestingly, the plateau in 
resistivity that was observed by Komiya et al}^ and cor- 
related with special behavior at 1/8 doping occurred at 
X = 0.13, which would correspond to p = 0.125 when 
corrected using our calibration curve. Note that, at all 
dopings, more than 85% of doped charge goes into the 
planar orbitals, though it appears increasingly in the api- 
cal bonds at higher doping. 



IV. CONCLUSIONS 

In summary, we have tested the hypothesis that dop- 
ing in La2-xSra;Cu04 is occurring predominantly in 
Cnd^z^^r'^-Opz orbitals^, rather than the canonical pic- 
ture of doping into the planar Cudx2_y2~Opx,y bands. 
Whilst the average structure evolves qualitatively as 
would be predicted according to Cnd^z^-r^-Opz dop- 
ing picture, more quantitative analysis suggests that the 
canonical picture of one hole doped into the planar bands 
per strontium is actually rather close to the true situa- 
tion. We present a correction factor quantifying the par- 
tition of doped charge between planar and apical bonds. 



6 



Acknowledgments 

We would like to acknowledge M. J. Gutmann, A. 
Bianconi, N. L. Saini and P. G. Radaelli for help and 
discussions. We acknowledge beamtime on the instru- 



ment GEM at ISIS. This work was supported by NSF 
through grant DMR-0304391. SJB acknowledges support 
and hospitality from Universita di Roma, "La Sapienza" , 
where part of the work was carried out. 



* Electronic address: |bozin@pa.msu.edu| 
URL: http://nirt.pa. msu Tedu/ 1 

1 J. K. Perry, J. Tahir-Kheli, and W. A. Goddard-III, Phys. 
Rev. B 65, 144501 (2002). 

2 D. Haskel, E. A. Stern, D. G. Hinks, A. W. Mitchell, and 
J. D. Jorgensen, Phys. Rev. B 56, R521 (1997). 

^ P. G. Radaelli, D. G. Hinks, A. W. Mitchell, B. A. Hunter, 
J. L. Wagner, B. Dabrowski, K. G. Vandervoort, H. K. 
Viswanathan, and J. D. Jorgensen, Phys. Rev. B 49, 4163 
(1994). 

* S. J. L. Billinge and M. G. Kanatzidis Chem. Commun. 
749 (2004). 

^ E. S. Bozin, S. J. L. Billinge, G. H. Kwei, and H. Takagi, 

Phys. Rev. B 59, 4445 (1999). 
® E. S. Bozin, Ph.D. thesis, Michigan State University, 2003. 

A.C. Larson and R.B. Von Dreele, Los Alamos Nat. Lab. 

Rep. LAUR 86-748 (1994). 

* T. Egami and S. J. L. Billinge, Underneath the Bragg 
Peaks: Structural analysis of complex materials, Pergamon 
Press, Elsevier, Oxford, England, 2003. 

^ P. F. Peterson, M. Gutmann, Th. Proffen, and S. J. L. 

Billinge, J. Appl. Crystallogr. 33, 1192 (2000). 
° Th. Proffen and S. J. L. BiUinge, J. Appl. Crystallogr. 32, 

572 (1999). 

^ S. Komiya, H.-D. Chen, S.-C. Zhang, and Y. Ando, Phys. 



Rev. Lett. 94, 207004 (2005). 

T. Hanaguri, C. Lupien, Y. Kohsaka, D.-H. Lee, M. 
Azuma, M. Takano, H. Takagi, and J. C. Davis, Nature 
430, 1001 (2004). 

R. D. Shannon and C. T. Prewitt, Acta Cryst. B 25, 925 
(1969). 

M. Braden, M. Meven, W. Reichardt, L. Pintschovius, 

M. T. Fernandez-Diaz, G. Heger, F. Nakamura, and T. Fu- 

jita, Phys. Rev. B 63, 140510 (2001). 

L D. Brown and D. Ahermatt, Acta Cryst. B 41 (1985) 

244. 

L D. Brown, J. Solid State Chem. 82, 122 (1989). 

J. Tahir-Kheli, (private communication). 

C. T. Chen, F. Sette, Y. Ma, M. S. Hybertsen, E. B. 

Stechel, W. M. C. Foulkes, M. Schluter, S-W. Cheong, 

A. S. Cooper, L. W. Rupp, Jr., B. Batlogg, Y. L. Soo, 
Z. H. Ming, A. Krol, and Y. H. Kao, Phys. Rev. Lett. 66, 
104 (1991). 

C. T. Chen, L. H. Tjeng, J. Kwo, H. L. Kao, P. Rudolf, 
F. Sette, and R. M. Fleming, Phys. Rev. Lett. 68, 2543 
(1992). 

B. Dalela, S. Dalela, N. L. Saini, R. K. Singhal, C. T. Chen, 
and K. B. Garg, Int. J. Mod. Phys. B 18, 2849 (2004). 



